Existence of the Abrikosov vortex state in two-dimensional type-II superconductors 

without pinning. 
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A nonperturbative fluctuation theory of the mixed state of two-dimensional type-II superconductor 
with finite size in the lowest Landau level approximation is proposed. The thermodynamic averages 
of the spatial average order parameter and of the Abrikosov parameter j3 a are calculated. It is shown 
that the thermodynamic fluctuations eliminate the Abrikosov vortex state in a wide region of the 
mixed state of two-dimensional type-II superconductor with real size. 



PACS numbers: 74.60. Ge 

Fluctuating behavior near second critical field, H C 2, 
in type-II superconductors has recently been under in- 
tense experimental |l|-|[ and theoretical |]6|-|l8|] study, 
particularly in connection with high temperature su- 
perconductors (HTSC). Investigation of the fluctuations 
changes step by step our habitual notion about nature 
of the mixed state. Concepts "vortex lattice melting" 
and "vortex liquid" have become very popular on 

first stage of a stormy investigation of fluctuation effects 
in HTSC. The vortex lattice melting transition is deter- 
mined experimentally by a change of resistive proper- 
ties in a perpendicular magnetic field This tran- 
sition was observed in bulk conventional superconduc- 
tors before the HTSC discovery |L9). It was interpreted 
in [i~9[ as a transition from the Abrikosov state into a 
one-dimensional state. But the vortex liquid as well as 
the one-dimensional state is not a new genuine thermo- 
dynamic state which is qualitatively different from the 
normal state ||. Therefore the interpretations by [jl]|| 
and by |ll| are just the same. The vortex lattice melt- 
ing is the transition from the Abrikosov state into the 
normal state with superconducting fluctuations (which 
was called as the one-dimensional state in |l^JT^]. The 
experimental investigations show that this transition oc- 



curs below H C 2 in all case [ 20| , |21| ,[1|-|3| . The second critical 
field line H C 2 (T) only marks a crossover from the normal 
metal to a strongly fluctuating superconducting state . 
I think that the widespread term " vortex lattice melting" 
is not enough right. Therefore I will use the term "tran- 
sition from the Abrikosov state into the normal state" in 
exchange for it. 

The starting point of vortex lattice melting theories is 
assumption of the existence of the Abrikosov state at low 
temperatures. This assumption bases oneself on experi- 
mental evidence only. But real samples with finite sizes 
are investigated experimentally whereas the thermody- 
namic limit is used and a influence of disorder (pinning 
centers) is not considered in the theory. Therefore the 
experimental evidence may be not enough. The pertur- 
bation fluctuation theory |2^,[l4| shows no sign of the 
Abrikosov state. Moreover, according to |ll| the fluctu- 
ation eliminates the Abrikosov state in superconductors 
with infinite sizes. The reduction of effective dimension- 



ality of the fluctuation in a magnetic field |2^] and the 
destruction of off-diagonal long range order by thermal 
fluctuations in the mixed state [^J^ JlT| , [r7|] cause doubt 
of the existence of the Abrikosov state in the thermody- 
namic limit also. According to results of p3] the fluc- 
tuation can not be considered as a perturbation in the 
mixed state of type-II superconductor with infinite sizes. 
Therefore a nonperturbative fluctuation theory of super- 
conductors with finite sizes is required. 

This theory for a most easy case is presented in this 
work. A two-dimensional superconductor in a perpendic- 
ular magnetic field near H C 2 is considered. In this case 
the lowest Landau level (LLL) approximation is valid and 
the system is zero-dimensional |lq]. Therefore the exact 
thermodynamic averages can be calculated easily. The 
theory is based on a new expression of the free energy 
obtained from the Ginsburg-Landau (GL) one, Fgl- The 
results obtained in this work are most close to one ob- 
tained in Q . In comparison with || a more easy method 
of calculation is suggested here and the concept of a tran- 
sition into the Abrikosov state of two-dimensional super- 
conductors with finite sizes is proposed. 

We will proceed from the accurate expression for the 
free energy 



F = -/c B Tln^exp( 



Fgl 



(1) 



Strongly type-II superconductors near H C 2 will be consid- 
ered there. In this case one can neglect the fluctuations 
in the magnetic field. The problem of fluctuations can 
be investigated within a framework of the GL free energy 
functional, with the order parameter confined to the low- 
est Landau level (LLL) for Cooper pairs. In this approxi- 
mation the GL free energy functional for two-dimensional 
superconductor may be write as (see Jf5|) 



F 



GL 



k B T 



y(e|*| 2 +0.5/3 a |*| 2 ) 



(2) 



e = Gi- 1 r\ht)- 1 / 2 {h-h c2 ), Gi = [2k B T/H*(0)(% 2 (0)} 
is the Ginsburg number in zero magnetic field; t = T/T c ; 
h = H/H c2 (0); H c2 (0) = -T c {dH c2 /dT) T = Te] H c (0) = 
—T c (dH c /dT); H c is the thermodynamic critical field; 
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£(0) is the coherence length at T = 0. d is the film thick- 
ness; V is the film area; |^| 2 = (J v d 2 r\ 1 i>\ 2 )/V is the spa- 

r 2 

tial average order parameter; (3 a — | X I/| 4 /| , J , | 2 is the gen- 
eralized Abrikosov parameter; l^l 4 = (J v d 2 r\^>\ A )/V . 
We use a dimensionless unit system in which l x = 1; 
l y /2 = 1; l x l y /2 = <P /H = 1; * = (dpUy^ksT) 1 / 4 ^. 
In this unit system the film area V is equal the degen- 
eracy number of the LLL. l x , l y are the parameter of the 
triangular Abrikosov lattice corresponding to a magnetic 
field value H; $o is the flux quantum; ip is order param- 
eter in the conventional units. To obtain (2) we used the 
well-known relation: a = —ehH C 2/mc; a 2 /2/3 = H 2 /8n. 
a and (3 are the conventional coefficients of GL free en- 
ergy 

In the LLL approximation the sum in (1) is to be taken 
over the subspace spanned by the LLL. To substitute the 
sum in (1) on two integral 



k B T) 

(la) 

N(\^\ 2 ; /3 a )c?|^I'| 2 c?/3 a is a subspace volume with the given 
values of |\&| 2 and (3 a . Pa is the minimum value of p a 
in the subspace spanned by LLL. It was shown in pi] ] 
that Pa correspond to the triangular lattice and is equal 
approximately 1.159595. 

To expand the order parameter on the Landau eigen- 
functions *(r) = V' 1 ' 2 ]T A f A (/> A (r). Then \W = 



F = -k B T\n[l d\% 



d/3 Q iV(|4'| 2 ;/3 a )exp(- 



V-'ExKR^x) 2 + (Im* A ) 2 ]. A 



eigenfunction number is equal the degeneracy number 
of the LLL, which is equal in our unit system V. 
Consequently the given values of \*f?\ 2 lie on a 2V- 

5 o 5 

dimensional sphere with radius \^\ 2 : (\^\ 2 ) 2 = 
J2x[( Re ^\/ v ° 5 ) 2 + (Im^x/V - 5 ) 2 ]. The area of this 
sphere is equal 27r y l/!(| , I'| 2 



T 0.5, 



'I 2 ) 



2V-1 



It is obvious that 



Pa value does not depend on the \^\ 2 value. Therefore 

iV(|*| 2 ;/3 a ) is proportional to |*| 2 °' 5 n(/3 Q ). n(fi a ) is 
a fraction of the 2V-dimensional sphere with the given 
value of Pa- Multipliers independent of |$| 2 and p a are 
omitted. 

To write the generalized Abrikosov parameter through 
eigenvalues 

A. = (Ea < ^ 1 ,a 2 ,A3,a 4 ^ 1 ^ 2 *A3*a 4 )/(EaI*a| 2 ) 2 . 
Va 1 ,a 2: a 3 ,a 4 = V' 1 J v d 2 r^l 1 (p* x ^x 3 (px 4 - We shall use 
the Eilenberger basis function [E7f. In this basis each 
ipx function describes the triangular Abrikosov lattice 
and A = (X x ;X y ) is their relative displacement. X x — 



A,, 



(l 2 y lL x )m. 



n; m are integer numbers; 



L x ;L y are film sizes. L x L y = V. Vx 1 ,x 2 .x 3 .^ 
lated in ||. 

Let *a=o = *q = 1. Following (24 
variables u+(A) and u_(A) 



was calcu- 



I introduce new 



(\Vo,o,x,-x \/2V ,o,x,-x) 1/2 lu+(\)+u_(\)} 



*-A = (\Vo^x,-x\/2V 0fl ^-x) 1/2 [u* + (X) -«* (A)] 

Then the generalized Abrikosov parameter may be 
write as 



Pa = Pa + 



£ A#0 (r + | u+ | 2 + r_M 2 ) + £/-i? 



[i + o.5£a #0 (KI 2 + M 2 )F 



(3a) 



T + = 2Vb,A,0,A + |Vb,o,A,-A| - Pa 



r_ = 2^0^,0^ - |Vo,o,a,~a | - Pa 



U =J2 ^ 1 ,A 2 ,A3,A 4 *A 1 *A 3 ^A 3 *> 
Ai#0 



B 



0a(£ I* 

A^O 
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The sum by all A ^ 0. 

At b « i, p a ~ A + EA^ ( r +l«+l 2 + r_M 2 ). 

Consequently at j3 a — Pa "C 1 the given value p a lie 
in a 2V-dimensional circle. A "length" of this circle is 
proportional to [(/?„ — Pa)°' 5 ] 2V ~ 2 - Therefore n(P a ) is 
proportional to (P a — Pa) v ~ x at P a — Pa <C 1. At a 
larger p a — Pa value n(P a ) increases more slowly with 
the P a — Pa increasing. In this case one may introduce a 
new function and write n(P a ) = [/(/3 a — Taking 
into account that V > 1 we may overwrite (la) as 



/>oo />oo 

F = -k B Tln[ / djl/j 2 / dp a exp 
Jo J/3 A 



where F^ ew is the new free energy expression in which 
the entropy term connected with |^| 2 and a values is 
take into account 

F Ne 



k R T 



= e|tf| 2 + 0.5/3 o |tf| 2 -ln|*| 2 -ln[/(/? a -/? A )] (4) 



Beca use the e xact t hermodynamic averages of the 

1 4" | 2 and p a values, < |\I/| 2 > and < p a >, are close to val- 
ues corresponded to the Fx ew minimum. Consequently 
these values are specified by relations 



<m 2 >=\( 



i 



11/2. 



2 < P a > < Pa> 



2<Pa> 



/'(< p a > -p A ) 
/(< Pa > -Pa) 



0.5 < I4-I 2 > 2 



(5a) 



(56) 
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The (5a) relation coincide with one obtained in |16| and 
is close to the entropy dependence obtained by O'Neill 
and Moore l|. This relation is close to the Abrikosov 
solution, < |^| 2 >= -e/ < p a >, at —e > 1 (below the 
H C 2 critical region) and to the result of the fluctuation 
theory in the linear approximation, < | \T/ 1 2 >= 1/e, at 
e 3> 1 (above the H C 2 critical region). 

The f(P a — Pa) function is not evaluated completely 
in this work. Therefore we calculate the thermodynamic 
averages of f3 a below the critical region only. At — e 3> 1, 
< P a > -j3 A < 1. Consequently /(< a > -f3 A ) ^< 
p a > —Pa and according to (5b) 



< P a > -P/ 



2 < m >- 2 ^2(^-) 2 

e 



20 A Giht(h c2 - h)~ 



Therefore in the cen- 
at t = 0.5, h = 0.25, 



The Ginsburg number values of conventional two- 
dimensional type-II superconductors lie in the interval 
3 * 10~ 6 < Gi < 3 * 10" 3 . 
ter of the mixed state region 
10- 5 < (< p a > -p A ) < KT 2 . 

According (lb) there is not thermodynamic singu- 
larity in two-dimensional type-II superconductor which 
may be connected with a phase transition from nor- 
mal state into Abrikosov vortex state. Therefore one 
should think that this transition occurs without thermo- 
dynamic singularity. I propose there a concept of this 
transition based on the difference between the Abrikosov 
state and the normal state with superconducting fluc- 
tuation in the Eilenberger function basis p7|| . It is 
obvious that above the H C 2 critical region all thermo- 
dynamic average of the Eilenberger eigenvalue squares, 

< I^aI 2 >, are equal and < ^(r)) 2 > is constant in 
space. Consequently, in the normal state with super- 
conducting fluctuation (in the vortex liquid phase j(|) 

< |$J2 >= y-i < |* A | 2 >=< |* A =o| 2 > whereas in 
the Abrikosov state < |*| 2 >= V^ 1 < |*a=o| 2 > in the 
mean field approximation. Therefore one should think 
that the transition into the Abrikosov vortex state take 
place when a one Eilenberger eigenvalue, \I>o > brings the 
main contribution to the order parameter: < IV'ol 2 >=< 

> /V < W > = (1 + Ea^o < IV'aI 2 >)- 1 - 



1. Where 



(1 + 0.5£a^o < \ u +\ 2 > + < \ u -\ 2 >)~ 
<|Va| 2 >=<|*a| 2 >/<|*o| 2 >. 

The < |-0o| 2 > value is connected with < (3 a > ~Pa 
value. We shall consider the region below the H C 2 critical 
region where < p a > ~Pa "C 1. Following to p4j one 
may show that T + ~ 20a and T_ ~ 3.1|A| 4 at |ApC 1. 
Consequently, according (3a), Ea#oI u +I 2 ^ 1 an d ^ 
may be that |it-| 2 S> 1 at j3 a — Pa <C 1. Therefore 
we will consider states with u + = only. In this case 
|V'o| 2 = (l + 0.5£ A#0 M 2 )- 1 and 



Pa ~Pa = 



A#0 



U-B 



(1 + 0.5£a^oM 2 ) 



(36) 



U = 0.25J2\^o Kxi,A 2 ,A 3 ,A 4 exp[0.5iOpAi + <^a 2 - Vx 3 ~ 
B = /3a(0.5£a^oM 2 ) 2 - Where exp(-i<p x ) = 

|Vb,0,A,-A|/Vb,0 : A,-A- 

At 0.5^ A ^ | w | 2 < 1 we have the Abrikosov state 

with fluctuation and at 0.5 ^2 



A#0 I 



| 2 > 1 we have 

the normal state with superconducting fluctuations. All 
states corresponding to very small p a — Pa value lie 
near the Abrikosov state (0.5 J2\=£o \ u - I 2 ^ ■*-)■ At 
larger f3 a — Pa value states exist which are removed from 
the Abrikosov state (states with 0.5^ A _^ |u_| 2 ^> 1). 
The critical value, (P a — Pa)o exists below which states 



(6) with 0.5 £ 



A#0 I 



" 3> 1 appear. I suppose that tran- 
sition into Abrikosov state occurs when < a > —Pa — 

{Pa-PA)c 

The {Pa — Pa)c value is not calculated in this paper. 
I will show only that in two-dimensional superconductor 
states with 0.5^ A _^ |it-| 2 ~S> 1 exist at small p a — Pa 
value. 

We will consider the states with U—(X X = 0;\ y ) = 
u/(X y L x ) = u/m at \X y \ < I < 0.1, m is odd number and 
U-(\ x ; X y ) = at all other A^, X y values. In this case 



Pa - Pa = ( 



0.52Pw 2 + (0.S7/Z - 2.8l)u 4 1 
(1 + 1.23u 2 ) 2 >T X 



(3c) 



At 0.5 Y,x^o \ u - 1 2 = L23u2 > !> Pa ~ Pa - (0.58// - 
1MI)/L X . The relation (3c) is valid at I < 0.1 only. At 
1 = 0.1, p a - Pa = 5.6/L x . Consequently (p a - p A ) c < 
5.6/L*. 

Thus the obtained results show that the Abrikosov vor- 
tex state can not exist in two-dimensional type-II super- 
conductors with infinite sizes (at L x — oo, (f3 a — Pa)c = 
and consequently < p a > —Pa > (P a — Pa)c always). 
The same result was obtained in P,fl6||. 

In "dirty" (Gi ~ 3 10~ 3 ) thin films with usual sizes 
(L x = L x /l x — 10 3 — 10 5 ) the transition into Abrikosov 
state takes place much below H C 2 (for example at t = 
0.5 and Gi = 3 * lO" 3 , < a > -p A < (p a - A ) C < 
5.Q/L X = 5*10" 3 -5*10" 5 below h = 0.25 ~ 0.5h c2 only, 
see above). It must be emphasize that these conclusion 
are valid for samples without pinning centers only. Any 
pinning of the flux lines by disorder, ets., will inhibit 
and slow the fluctuations and will stabilize the Abrikosov 
state |J. 

According to the concept described above, the 
Abrikosov state and the normal state with supercon- 
ducting fluctuation have different distributions of the or- 
der parameter in the space. Therefore a transition from 
the Abrikosov state into the normal state (vortex lat- 
tice melting) can be observed experimentally by pinning 



3 



appearance (by a change of resistive properties in a per- 
pendicular magnetic field). It was observed in that 
this transition occurs (pinning appears) much below H C 2 
in thin amorphous films with a small amount of pinning 
centers. It may be considered an experimental confirma- 
tion of the result obtained above. 
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